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Abstract. We analysed data from five XMM-Newton observations of GX13+1 to investigate the variability of 
the photo-ionised absorber present in this source. We fitted EPIC and RGS spectra obtained from the "least- 
variable" intervals with a model consisting of disc-blackbody and blackbody components together with a Gaussian 
emission feature at ~6.55-6.7 keV modified by absorption due to cold and photo-ionised material. We found a 
significant correlation between the hard, ~ 6-10 keV, fiux, the ionisation and column density of the absorber 
and the equivalent width of the broad iron line. We interpret the correlation in a scenario in which a disc wind 
is thermally driven at large, ~ 10^" cm, radii and the broad line results from reprocessed emission in the wind 
and/or hot atmosphere. The breadth of the emission line is naturally explained by a combination of scattering, 
recombination and fluorescence processes. We attribute the variations in the absorption and emission along the 
orbital period to the view of different parts of the wind, possibly located at slightly different inclination angles. 
We constrain the inclination of GX 13-f 1 to be between 60 and 80° from the presence of strong absorption in the 
line of sight, that obscures up to 80% of the total emission in one observation, and the absence of eclipses. We 
conclude that the presence of a disc wind and/or a hot atmosphere can explain the current observations of narrow 
absorption and broad iron emission features in neutron star low mass X-ray binaries as a class. 

Key words. X-rays: binaries - Accretion, accretion disks - X-rays: individual: GX 13-1-1 
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1. Introduction 



GX 13-1-1 is a type I X-ray burster (|FleischmaDl 11985 ; 
iMatsuba et al. 1995 ) with a bright persistent X-ray emis- 
sion. It is associated with an evolve d late-type K5 iii star, 
locat ed at a distance of 7±1 kpc ( Bandvopadhvav et al.l 
Il999l) . GX13-f-l is the most luminous "atoll" source 
known to date, with an intermediate luminosity ('-~^0.1 
Eddington luminosity) between the so-called "Z" Low 
Mass X-ray Binaries (LMXBs) and all other "atoll" 
sources. Its nature has been debated in the literature 
since it s hows some features from Z sources, like 57- 
69 QPOs (iHoman et al.l 1998 ). persistent radio emission 
(|Grindlay fc Seaauistlll986l ) and high accretion rate, and 
some proper from atoll sources (e.g. the path in the colour- 
colour diagram). The investigation of the ape riodic vari 



ability of bright persistent X-ray sources led iReig et al 



([20031) to propose a new classification for LMXBs, depend- 
ing on the rms amplitude and the slope of their power den- 
sity spectra. In their scheme, GX 13+1 is the only "atoll" 
source which has been classified as a type I source (lower 
rms, rms<20%, and flatter power spectra, a < 0.9), a class 



which includes all the "Z" sources. ISchnerr et al.l (|2003f ) 



proposed to explain the source peculiar spectral and power 
spectral properties with the presence of a relativistic jet. 
[Paizis et al. , (2006.) first detected a hard X-ray tail in the 
GX 13-1-1 average hard X-ray spectrum (>20 keV), sum- 
ming together all the available INTEGRAL observations. 
They found a correlation for "Z" sources and bright "GX" 
atolls (like GX13+1), between the X-ray hard tail (40- 
100 keV) and their radio luminosity, likely indicating that 
hard tails and energetic electrons causing the radio emis- 
sion may have the same origin, most likely the Compton 
cloud located inside the neutron star (NS) magnetosphere. 



Corbetl (|2003f) found a modulation of 24.065±0.018 d 



in RXTE ASM data spanning 7 years and suggested that 
this may be the orbital pe riod of the system. A s imilar 
modulation was found by iBandvopadhvav et al. ( 2002 ) 
after performing infrared photometry of the source. The 
24 day period has been recently confirmed f rom both the 
K-ba nd and additional X-ray o bservations (ICorbet et aJ " 



2010l ). Other NIR observations (jProning et al.ll2008h show 



that the accretion disc is the dominant emitter in GX 13-1-1 
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(~3/4 of NIR flux comes from the disc at the time of those 
observations) . 



2. XMM-Newton observations 



Ueda et aLl (|2001l ) discovered narrow resonant ab- 
sorption hnes near 7 keV from GX13+1 using ASCA. 
The same f e ature s were studied with XMM-Newton by 
Sidoh et all (|2002l ) who found a complex of features iden- 
tified with resonant scattering from the Ka and K/3 lines of 
He- and H-hke iron (FeXXV and Fexxvi), H-hke calcium 
(Caxx) Ka as well as the presence of a deep (r ~ 0.2) 
Fe XXV absorption edge at 8.83 keV. The absorption lines 
were superposed on a broad emission feature whose en- 
ergy and breadth were poorly determined, partly due to 
the presence of the deep Fe Kg features. Chandra HETGS 
observations ( Ueda et al.ll2004l ) revealed that the features 
are blue-shifted indicating an outflowing plasma with a 
velocity of ~400 km s~^. 

GX13+1 is unusual in that all the other LMXB sys- 
tems that exhibit prominent FeXXV and FeXXVi feature s 
are dipping sources (see Table 5 of iBoirin et al.l 120041 ). 



Dipping sources are LMXB systems that are observed 
from close to the plane of the accretion. The dips are be- 
lieved to be caused by periodic obscuration of the central 
X-ray s ource by structure loc ated in the outer regions of 
a disc (jWhite fc Swank! 1 19821 ). The depth, duration and 
spectral evolution of the dips varies from source to source 
and often from cycle to cycle. The lack of an orbital mod- 
ulation of the absorption features outside of the dips in- 
dicates that the absorbing plasma is located in a cylindri- 
cal geometry around the compact object. The important 
role that photo- ionised plasm as pl ay in LMXBs was reco; 
nised by [SoirineEaD (12005*) and' Pi'az Trigo et al.l (|200 
who were able to model the changes in both the narrow 
X-ray absorption features and the continuum during the 
dips from all the bright dipping LMXB observed by XMM- 
Newton by an increase in the column density and a de- 
crease in the amount of ionisation of a photo-ionised ab- 
sorbing plasma. Since dipping sources are normal LMXBs 
viewed from close to the orbital plane, this implies that 
photo-ionised plasmas are common features of LMXBs. 
Outside of the dips, the properties of the absorbers do 
not vary strongly with orbital phase suggesting that the 
ionised plasma has a cylindrical geometry with a maxi- 
mum column density close to the plane of the accretion 
disc. 

Since GX 13-1-1 exhibits deep Fe absorption features 
similar to the dipping LMXBs, this suggests that GX 13-t-l 
may be a high-inclination LMXB and that as such the 
overall X-ray continuum may be similarly affected. To in- 
vestigate the origin of the variability of the absorption and 
the presence of dipping activity, we successfully applied for 
five observations of GX 13-1-1 with XMM-Newton, spaced 
by at least one day. In this paper, we report on this mon- 
itoring program. We include results from the RGS detec- 
tors and the optical monitor, which were not reported for 
the early 2000 observations. 



The XMM-Newton Observatory (jJansen et al.l 120011) in- 
cludes three 1500 cm^ X-ray telescopes each with an 
EPIC (0.1-15 keV) at the focus. Tw o of the EPIC imag- 



ing spectrometers use MOS CCDs (iTurner et al.l 120011) 



and one uses p n CCDs (IStriider et al.l 20011) . The RGSs 



(0.35-2.5 keV, iDen Herder et al.l l200lD are located be- 
hind two of the telescopes. In addition, there is a co- 
aligne d 30 cm diameter Optical/UV Monitor telescope 
fOM. iMason et al.l[200l . providing simultaneously cov- 
erage with the X-ray instruments. Data products were 
reduced using the Science Analysis Software (SAS) ver- 
sion 10.0. The EPIC MOS cameras were not used during 
the observation in order to allocate their telemetry to the 
EPIC pn camera and avoid Full Scientific Buffer in the lat- 
ter. We present here the analysis of EPIC pn data, RGS 
data from both gratings and OM data. 

Table |T] is a summary of the XMM-Newton observa- 
tions. We used the EPIC pn in timing mode. In this mode 
only one CCD chip is operated and the data are collapsed 
into a one-dimensional row (4.'4) and read out at high 
speed, the second dimension being replaced by timing in- 
formation. This allows a time resolution of 30 /is. We used 
the SAS task epf ast on the event files to correct for a 
Charge Transfer Inefficiency (CTI) effect which has been 
observed in EPIC pn timing mode when high count rates 
are presenfl Ancillary response files were generated us- 
ing the SAS task arf gen following the recommendations 
of the XMM-Newton SAS User guide for piled-up obser- 
vations in timing mode, whenever applicable. Response 
matrices were generated using the SAS task rmf gen. 

Light curves were generated with the SAS task 
epiclccorr, which corrects for a number of effects like 
vignetting, bad pixels, PSF variation and quantum effi- 
ciency and accounts for time dependent corrections within 
a exposure, like dead time and GTIs. 

The SAS task rgsproc was used to produce cali- 
brated RGS event lists, spectra, and response matrices. 
We also chose the option keepcool=no to discard single 
columns that give signals a few percent below the values 
expected from their immediate neighbours. Such columns 
are likely to be important when studying weak absorption 
features in spectra with high statistics. We used the SAS 
task rgsbkgmodel to compute model background spectra 
from RGS background templates. We generated RGS light 
curves with the SAS task rgslccorr. 

The OM was operated in image+fast mode. The U, 
UVWl and UVM2 filters were used. In this mode the in- 
strument produces images of the entire 17' x 17' FOV 
with a time resolution between 800 and 5000 s and event 
lists with a time resolution of 0.5 s from a selected 11" x 
11" region. The SAS task omf chain was used to extract 



^ More information about the CTI correction can be found 
in the EPIC status of calibration and data analysis and in the 
Current Calibration File (CCF) release note Rate- dependent 
CTI correction for EPIC-pn timing modes, by Guainazzi et al. 
(2009), at http://xmm.esac.esa.int/external/xmm_calibration 
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Table 1. XMM-Newton observations of GX 13+1. Results from the first three observations are reported in lSidoh et al 



( 20021 ). T is the total EPIC pn exposure time and C is the pn 0.6-10 keV total count rate after dead time correction. 
In all cases the EPIC pn thin filter was used. For obs 5 only RGS and OM data are available. 



Obs 


Observation 


Observation Times (UTC) 


T 


C 


Num 


ID 


.Start 




End 


(ks) 


is-') 






(year month day) 


(hr:mn) 


(hr:mn) 






1 


0122340101 


2000 March 30 


14:10 


16:58 


6.0 


674 


2 


0122340901 


2000 April 01 


04:29 


07:54 


8.1 


781 


3 


0122341001 


2000 April 01 


08:51 


11:16 


4.6 


690 


4 


0505480101 


2008 March 09 


15:48 


22:50 


22.9 


843 


5 


0505480701 


2008 March 11 


19:00 


23:10 


12.0 




6 


0505480201 


2008 March 11 


23:10 


27:16 


12.4 


909 


7 


0505480301 


2008 March 22 


02:20 


07:15 


15.3 


818 


8 


0505480501 


2008 March 25 


23:01 


26:59 


11.9 


761 


9 


0505480401 


2008 September 5 


21:49 


26:48 


15.5 


746 
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Fig. 1. a) 0.6-10 keV EPIC pn (dots) and 0.6-1.8 keV RGS (squares) background subtracted light curves for each 
GX13+1 observation with a binning of 64 s. The left (right) axes show the scale for EPIC pn (RGS). b) Hardness 
ratio (counts in the 2-10 keV band divided by those between 0.6-2 keV) for the EPIC pn light curves. 



light curves of GX 13+1 from the high time resolution fast 
mode data. We used a time sampling of 100 s in the curve 
extraction to improve the signal-to-noise ratio. 



2.1. Pile-up and X-ray loading in the EPIC pn camera 

The count rate in the EPIC pn was close to, or above, the 
800 counts s~^ level, at which X-ray loading and pile-up 
effects become significant. 
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Pile-up occurs when more than one photon is read in 
a pixel during a read-out cycle. This causes photon loss, 
pattern migration from lower to higher pattern types and 
hardening of the spectrum, because the charge deposited 
by more than one photon is added up before being read 
oulH. In addition, when high count rates are present, the 
offset map calculated at the beginning of a exposure may 
be contaminated by X-ray events from the source, the so- 
called "X-ray loading" . As a consequence pattern migra- 
tion from higher to lower pattern types and a shift to 
lower energy for all the events associated with the contam- 
inated pixel occujfl. Since both X-ray loading and pile-up 
cause significant spectral distortion, we investigated in de- 
tail their presence before extracting the spectra. 

We used the SAS task epatplot, which utilizes the 
relative ratios of single- and double-pixel events which de- 
viate from standard values in case of significant pile-up, as 
a diagnostic tool in the pn camera timing mode data and 
found that all the spectra were affected by pile-up. Next, 
we extracted several spectra selecting single and double 
timing mode events (patterns to 4) but different spa- 
tial regions for the source. Source events were first ex- 
tracted from a 62" (15 columns) wide box centred on the 
source position (Region 1). Next we excluded 1, 3, 5 and 
7 columns from the centre of Region 1 (Regions 2-5) and 
extracted one spectrum for each of the defined regions. 
We obtained spectra free of pile-up once events within the 
inner 5 columns were excluded. 

Next, we inspected the pn offset maps to diagnose the 
presence of X-ray loading. For each observation, we cal- 
culated the so-called "residual" offset map by subtracting 
the offset map of a nearby observation taken with closed 
filter from the offset map of the observation. We found 
that all the observations showed effects of X-ray loading 
in at least the inner 5 columns of the PSF. X-ray loading 
could be significant in the inner 7 columns, but the noise 
fluctuations in the individual rows for outer columns are of 
the order of the level of excess in the 6th and 7th columns. 
Therefore, we extracted the final spectra after excluding 
events from the central 5 columns. 



2.2. Background subtraction 

In the EPIC pn timing mode, there are no source-free 
background regions, since the PSF of the telescope extends 
further than the central CCD boundaries. The central 
CCD has a field of view of 13.'6 x 4.'4 in the pn. In timing 
mode, the largest column is the one in which the data are 
collapsed into one-dimensional row. Therefore, the maxi- 
mum angle for background extraction is 2', compared to 
5' for imaging modes. Since GX 13+1 is very bright, its 
spectrum will not be significantly modified by the "real" 



background which contributes less than 1% to the total 
count rate in most of the bandwidth. Conversely, subtract- 
ing the background extracted from the outer columns of 
the central CCD will modify the source spectrum, since 
the PSF is energy dependent and the source photons scat- 
tered to the outer columns do not show the same energy 
dependence as the photons focused on the inner columns. 
Therefore, we chose not to subtract the "background" ex- 
trac ted from the outer regio i is of the centra l CCD (see 
also lDone fc Di'az TrigollioioUNg et al.ll2010l ). We expect 
a contribution from the background to the total count rate 
of more than 1% below 1.7 keV. Therefore, below this en- 
ergy we used only the RGS spectra, for which background 
templates were available. 



3. Light curves 

3.1. EPIC pn and RGS light curves 

Figurefflshows 0.6-10 keV EPIC pn and 0.6-1.8 keV RGS 
light curves of five XMM- Newton observations of GX 13-1-1 
performed in 2008 together with three observations per- 
formed in 2000, for comparison, with a binning of 64 s. 

The EPIC pn light curves are remarkably similar to 
the simultaneous RGS light curves despite the different 
energy band. We note that the lower RGS count rate with 
respect to the pn count rate in the 2008 observations com- 
pared to the 2000 observations is due to the inoperative 
CCDs 7 (RGSl) and 4 (RGS2) since early after the be- 
ginning of the mission. A high variability is present in all 
the observations. The average count rate changes signifi- 
cantly from one observation to the next one, with values 
between ^^700 and 900 counts s~^. The count rate vari- 
ability is not associated to significant changes in the spec- 
tral hardness for obs 4 and 6. However, obs 1, 3, 8 and 
9 show spectral hardening associated with a decrease of 
count rate. Obs 6 shows the highest count rate, with peaks 
of up to ~1100 counts s^^ in the EPIC pn camera. Obs 8 
shows the lowest count rate, with less than 1 counts s^^ 
at the beginning of the observation in the RGS, compared 
to ~6 counts at the end of the observation. This be- 
haviour could be related to dipping activity. Obs 9 shows 
a decrease of count rate at ~23-24.5 h which is related to 
spectral hardening and could be again an indication for 
dipping behaviour in GX 13-1-1. 

We next examined the hardness ratio (counts in the 
2-10 keV band divided by those between 0.6-2 keV) as a 
function of 0.6-10 keV count rate for all the observations 
(Fig. [5]). The largest hardness ratios correspond to the 
intervals of low count rate present in obs 1, 3, 8 and 9. This 
behaviour is typical of deep dipping states of dippers. A 



See XMM-Newton Users Handbook for more information 
on pile-up. 

^ More information about X-ray loading can be found in the 
PN X-ray loading investigation results calibration document by 
Smith (2004) at http://xmm.esac.esa.int 



detailed comparison of Fig.[2]with Fig. 4 from lBoirin et al 
(|2005l ) shows that indeed obs 1, 3, 8 and 9 occupy the same 



parameter space in the figure as the deep dipping state of 
4U 1323— 62, while obs 2, 4 and 7 resemble the shallow 
dipping states and obs 6 the persistent state. We analyse 
in detail the origin of these states in Sect. 14.1.41 
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Fig. 2. Hardness ratio (2-10 keV/0.6-2 keV counts) ver- 
sus 0.6-10 keV EPIC pn count rate for all the observations. 
Data from obs 1 to 3 are shown as filled squares, filled tri- 
angles and filled circles, respectively. Data from obs 4 and 
6 to 9 are shown as open squares, open triangles, open cir- 
cles, crosses and crooked crosses, respectively. Each point 
corresponds to a binning of 64 s. 



3.2. Light curves from other X-ray missions 

In order to investigate the variability of GX 13+1 further, 
we extracted from the High Energy Astrophysics Science 
Archive Research Center (HE AS ARC) light curves from 
previous GX 13-t-l X-ray observations with the maximum 
time resolution available. We did not perform a systematic 
search of dips in all the GX13+1 archival observations. 
Instead, we looked for the signature of dipping episodes in 
observations where a "bi-modal" behaviour of the count 
rate with hardness ratio had been reported in the litera- 



ture ( Stella et al. 1985: Schnerr et al 



20031 ). Our goal was 



to determine if dipping behaviour may have been not rec- 
ognized in previous observations of GX 13+1 and to see 
if such dipping behaviour, if existent, is associated with a 
certain state of the source. 

Details on the observations of GX13+1 from archival 
X-ray observations which may show dipping behaviour are 
given in Table [5] The light curve corresponding to the 
RXTE observation is shown in Fig. [3] The light curve 
from the EXOSA T arch ival observation is shown in Fig. 1 
(Il985h . 



from Stella et al 



In both observations a decrease 
of count rate, which is more pronounced in the soft band, 
is visible. The decrease in the 1-4 keV EXOSAT (2-4 keV 
RXTE) light curves is ~25-30% and is accompanied by 
spectral hardening. 

Finally we examined the GX13+1 RXTE All-Sky 
Monitor (ASM) 1.5-12 keV light curves spanning more 
than 13 years and with a time resolution of 96 s to search 
for any characteristic signature at the time of appearance 
of dips. We corrected the ASM light curves to the solar sys- 
tem barycentre using the FTOOL task f axbary and folded 
the light curves into segments between 0.1 and 24 days 
duration, but clear dipping activity was not detected. A 




11.4 11.5 11.6 11.7 

Time (hours ol 1998 May 17) 



Fig. 3. 2-4 keV (middle) and 4-9 keV (bottom) RXTE 
PCA background subtracted light curves for an archival 
GX13+1 observation with a binning of 16 s. Top: 
Hardness ratio (counts in the 4-9 keV band divided by 
those between 2-4 keV) for the RXTE PCA light curves 
shown in the middle and bottom panels. 

Table 2. X-ray archival observations of GX 13+1 with 
potential dip states. 



Date 


Mission 


Reference 


1983 September 22 


EXOSAT 


Stella et al. (19851 


1998 May 17 


RXTE 


Schnerr et al. f2003) 



modulation of the ASM data has been alrea dy extensively 
searched and it is r eported in the literature (|Corbetil2003 : 
Corbet et al1l2010l ). Therefore, we did not look further for 
periodic behaviour in the ASM data. 



4. X-ray spectra 

We extracted EPIC pn and RGS spectra for each observa- 
tion. We rebinned the EPIC pn spectra to over-sample the 
FWHM of the energy resolution by a factor 3 and to have 
a minimum of 25 counts per bin, to allow the use of the 
statistic. To account for systematic effects we added 
a 0.8% uncertainty to each spectral bin after rebinning 
(we note that this is equivalent to add ^2% uncertainty 
to each spectral bin before rebinning). We used the RGS 
spectra with two different binnings: in Sect. 14.11 we re- 
binned the RGS spectra to have a minimum of 25 counts 
per bin, to be able to consistently use the statistic 
for both EPIC pn and RGS spectra. In Sect, we re- 
binned the RGS spectra to over-sample the FWHM of 
the energy resolution by a factor 3 to be sensitive to nar- 
row features and we used the C-statistic (Cash' '1979"). We 
performed spectral analysis using XSPEC (Arnaud,,1996[l 
version 12.6.0q. Since there were no photons in the RGS 
spectra below ^0.8 keV, we used the RGS spectra in the 
energy interval 0.8-1.8 keV. We used the pn spectra be- 



6 



Di'az Trigo et al.: GX13+1 



tween 1.7 and 10 keV in order to exclude energy bins for 
which we expect the spectrum to be affected by back- 
ground (see Sect. 12. 2p . To account for absorption by neu- 
tral gas we used the tbabs XSPEC model with solar abun- 
dances ( Anders fc Grevess3ll989l) . Spectral uncertainties 
are given at 90% confidence (Ax^= 2.71 for one interest- 
ing parameter), and upper limits at 95% confidence. 

We extracted intervals of low variability for all the ob- 
servations, using as a criterion the hardness ratio and the 
count rate of the observations. For obs 7 there are no sig- 
nificant variations of hardness ratio or count rate during 
the observation. Therefore we only consider one interval 
for this observation. We divided obs 4, 6, 8 and 9 in two 
intervals each ( "high" and "low" ) and examined the spec- 
tra separately. A quick look shows that while for obs 4, 
8 and 9 both intervals show significant spectral changes, 
the two intervals of obs 6 are only different in the normal- 
isation factors. Therefore we continued the analysis for 
obs 6 only with one interval. We continued the analysis 
for obs 4, 8 and 9 only with the "high" , less variable, in- 
tervals, for which we could accumulate a sufficient, > 5 ks, 
net exposure time. 

4.1. EPIC and RGS spectral analysis 
4.1.1. Phenomenological model 

We first examined the spectra of the low-variability in- 
terval separately for each observation. We combined the 
RGS spectra (RGSl and RGS2, order 1) and the EPIC 
pn spectrum and fitted them with a model consisting of 
a disc blackbody and a blackbody components modified 
by neutral material. A constant factor, fixed to 1 for the 
pn spectrum, but allowed to vary for each RGS spectrum, 
was included multiplicatively in order to account for cross- 
calibration uncertainties. 

The residuals showed strong absorption lines near 
7 keV superposed on a broad emission feature and weak 
absorption features at ^1.8 and 2.3 keV. The latter are 
most likely due to residual calibration uncertainties (after 
using the SAS task epf ast), due to an incorrect applica- 
tion of the CTI correction in the EPIC pn camera when 
high count rates are presenl|f|. We modeled such features 
with two Gaussian absorption components and do not dis- 
cuss them further. However, we note that these features, 
if due to residual uncertainties of the CTI correction, in- 
dicate that the energy gain could be compromised in the 
whole energy band. Detailed plots of the residuals from 
the best-fit continuum model around the Fe K region are 
shown in Fig. 01 Absorption features from highly ionised 



* We note that while the use of epf ast is recommended 
by the EPIC calibration team, the recent discovery that the 
count-rate dependency of the CTI may be caused by X-ray 
loading implies that the epf ast correction may improve, but 
not completely solve, the CTI calibration deficiency (for more 
details, see "Evaluation of the spectral calibration accuracy 
in EPIC-pn fast modes" at http://xmm.esac.esa.int/external/ 
xmm_calibration) . 



species of iron, such as Fexxv and Fexxvi Ka and K/3 
are evident. A broad iron emission line is present in all 
the observations. Both the absorption and emission fea- 
tures show significant variations among observations, with 
obs 6 and 7 showing the most ionised species in absorp- 
tion and the weakest line in emission. Obs 4 and 9 show 
the strongest line in emission and a prominent absorption 
edge at ~9 keV, that we identify with absorption from 
Fexxv and Fexxvi. 

We first evaluated the effect of the continuum in the 
absorption and emission features in a mo del- independent 
manner calculating the ratio among the spectra for all 
the observations with respect to obs 6 (see Fig. [5]). Obs 7 
shows a similar spectral shape to obs 6 and a flux lower by 
~10% in aU the energy band. Above ~6 keV, obs 4, 8 and 

9 show a significant deficit of photons with respect to obs 6 
and 7. Below ~5 keV, obs 8 and 9 show the lowest fluxes. 
Interestingly, obs 4 shows the lowest 6-10 keV flux of all 
the observations but the highest 1.5-5 keV flux together 
with obs 6. This analysis indicates that the hard, 6- 

10 keV, flux could be correlated with the ionisation of 
the absorption features and the significance of the broad 
emission line. As the hard flux increases from obs 4 to 
obs 9, 7 and 6 the significance of the broad emission line 
decreases and the absorption features become weaker and 
reveal an increase of ionisation in the plasma. Obs 8 shows 
an intermediate hard fiux and broad emission line, but the 
less ionised absorption features, indicating that the flux 
below ~5 keV may play an important role in the ionisation 
of the absorbing plasma. 



0.6 




Fig. 5. Ratio of EPIC pn GX 13+1 spectra from obs 4, 7, 
8 and 9 with respect to obs 6. 

In order to quantify the relation between the sig- 
niflcance of the emission and absorption features and 
the changes in the continuum, we next included in 
the model Gaussian absorption and emission features 
and absorption edges to account for the residuals near 
7 keV. The total model consisted of disc blackbody and 
blackbody components, one Gaussian emission feature at 
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Energy [keV] Energy [keV] 

Fig. 4. Ratio of the data to the continuum model at the Fe K region for obs 4 (upper-left), 6 (upper-middle), 7 
(upper-right), 8 (lower- left) and 9 (lower-middle). The dotted vertical lines indicate from left to right the rest energy 
of the transitions of neutral Fe, Fexxv Ka, Fexxvi Ka, Nixxvii , Fexxv K/3, Nixxviii and Fexxvi K/3. 



~6.6 keV and four Gaussian absorption features at the 
energies of Fexxv and Fexxvi Ka and Kf3, modified 
by photo-electric absorption from neutral material and 
by absorption edges at the energies of Fe xxv and Fe xxvi 
(tbabs*edgei*edge2* (diskbb+bbodyrad+gaui+gau2+gau3 
+gau4+gau5), hereafter Model 1). The fits with Model 1 
were acceptable, with xl between 1.0 and 1.2 for ~530- 
700 d.o.f. for obs 4-9. The parameters of the Gaussian 
emission and absorption features and edges for the best-fit 
with this model are given in Table [3] 

We can infer some properties of the absorber from the 
equivalent width (EW) of the absorption features, the 
depth of the edges and the predominance of a given ion. 
For example, the depth of the Fe xxv edge indicates that 
obs 4 and 9 have a column density of ionised plasma sig- 
nificantly larger than obs 6-8. The ratio of the Fexxv 
and Fe xxvi Ka to K/3 lines in obs 4 shows saturation of 
the lines, indicating most likely a column density above 
10^^ atom cm~^ 0. The ratio of the EW of the Fexxv 
to Fexxvi lines determines the ionisation of the plasma 
and shows that obs 6 and 7 have a very ionised plasma 
and that obs 8 has the least ionised plasma. Confirmation 
of the ionisation is given by the presence of low signifi- 
cance absorption features at ~1.88 keV in obs 8 and 9, 
that we identify with Sixiii, and at ^2.62 keV in obs 4, 8 
and 9, that we identify with Sxvi. Obs 8 shows also small 
features at ~2.0, 2.5 and 4.2 keV that we identify with 
Si XIV, Sxv and Caxx. Taking the above into account we 
can also examine the relation between emission and ab- 

^ The exact column density at which the absorption lines 
reach saturation depen ds on the velocity broadening of the 
lines (see e.g. Fig. 1 of iTombesi et "alll201ll ). 



sorption lines. The EW of the broad emission line is cor- 
related with the EW of the narrow absorption lines (see 
Fig. in]). Obs 6 and 7 show the smallest emission line and 
the thinnest absorbing plasma. Obs 4 and 9 show instead 
the largest emission lines together with the largest col- 
umn density of ionised plasma. The EW of the emission 
and absorption lines for obs 8 are intermediate between 
obs 4 and 9 on one side and obs 6 and 7 on the other 
side. Finally, the column density of the ionised plasma in- 
creases as the plasma becomes less ionised for obs 4, 6, 
7 and 9. This is expected since obs 6 and 7 are likely to 
have a fraction of fully ionised plasma that reveals itself 
only through electron scattering. Instead, for obs 8 the 
decrease in the column density with respect to obs 4 and 
9 cannot be explained simply by a change in ionisation 
and indicates that at the time of obs 8 a smaller column 
density is present. 

4.1.2. Photoionised plasma model 

The strength of the absorption features shown in Table [3] 
indicates that the continuum may be significantly affected 
by the absorbing plasma. Therefore, to quantify the above 
changes in a more physical manner, we substituted the 
absorption features and edges by the component warmabs, 
which models the absorption due to a photoionised plasma 
in the line of sight. This component not only accounts 
for the narrow absorption features evident near 7 keV 
but also modifies the overall continuum shape at regions 
where the spectral resolution is not enough to resolve in- 
dividual features. We note that the warmabs model does 
not include Compton scattering, in contrast to the xabs 
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Table 3. Parameters of the most significant absorption lines and edges and of the broad emission line detected in 
the EPIC pn spectra of obs 4-9. / indicates that a parameter was fixed. The width of the absorption features was 
constrained to be ^0.1 keV and of the broad emission feature to ^ 1 keV. p indicates that the error of a parameter 
pegged at the limit imposed, gaui corresponds to the broad emission feature and gau2-gau5 to narrow absorption 
features. The energy of the FeXXV feature of obs 7 was unconstrained due to its low significance. Therefore, we fixed 
its energy considering that it would show a similar blueshift as the Fe xxvi feature for the same observation. 



Observation No. 



Parameter 



Comp. 



£;gau (keV) (i-53tom 6.63 ±0.13 6-60 toTi 6.65 1:^ ?^ 6.71 ±0.08 

(7 (keV) 0.72to.lt 0.74±0.19 0.62tol9 0.61±o'.18 O-GOto.H 

fcgau (10"^ ph cm"^ s~^) 12 ±4 8±2 6±3 9±2 11 ±2 

EW (eV) 185 ± 62 108 ±27 84 ± 42 139 ± 31 177 ± 32 
gau2 

£gau (keV) Q-7^to°ol - 6.73(f) 6.71 ±0.02 6.73 ±0.02 

a (keV) < 0.04 - 0.008 laoos" < 0-07 < 0.05 

fcgau (10"" ph cm"^ s"^) 12±3 - 4±3 16tt U±i 

EW (eV) 20±5 - 6±5 241^2 21±6 
gau3 

£gau (keV) 7.02 ±0.06 7.03 ±0.02 7.00 ±0.01 7.02 ±0.02 7.03 ±0.01 

CT (keV) <0.03 0.07 ±0.05 < 0.06 <0.04 <0.06 

fegau (10"" ph cm"^ s"^) 23 ±4 20 ± 5 22 ±3 18 ±4 24 

EW (eV) 41 ± 7 30 ± 8 34 ± 5 32 ± 7 42 tli 
gau4 

Sgau (keV) 7.9lt°oil - - 7.91 ±0.05 7.84+°Z 

a(keV) 0.003 iHo? " " 0.0051^^0? 0.03 

fcgau (10"* ph cm"^ s"i) 11 ±4 - - 5±3 5±4 

EW (eV) 33 ± 12 - - 15 ± 9 13 ± 10 

Sgau (keV) 8.24:toot " 8.20±0.10 8.34±0.05 8.26±0.04 

a(keV) 0.1 - 0.1 1^, < 0.097 0.02 

fcgau (10"* ph cm~^ s"^) 10±4 - 6±3 5±2 7±3 

EW (eV) 37±15 - 15±8 18 ± 7 25±11 
edgei 

E (keV) 8.83 (f) 8.83 (f) 8.83 (f) 8.83 (f) 8.83 (f) 

r 0.18 ±0.03 <0.04 <0.04 <0.05 0.08 ±0.04 
edge2 

E (keV) 9.28 (f) 9.28 (f) 9.28 (f) 9.28 (f) 9.28 (f) 

r <0.03 0.04 ±0.03 0.04 ±0.03 < 0.05 <0.09 

xlid-o.i.) 1.08 (574) 1.11 (694) 1.03 (678) 1.07 (530) 1.18 (561) 



model in the SPEX package (jKaastra et al.Ml996f ). The 
contribution of Compton scattering is expected to be sig- 
nificant when modeli ng the changes between persistent 



and d ipping emission (iBoirin et al.ll2005t iDiaz Trigo et al 
20061) . Therefore we added the component cabs to the 



model, which accounts for non-relativistic, optically-thin 
Compton scattering. We forced the column density of 
the cabs model to be equal to the column density of 
the warmabs component multiplied by a factor of 1.21. 
The latter factor accounts for the number of electrons 
per hydrogen atom for a material of solar abundances 



(|Stelzer et all Il999l) . The non-relativistic approximation 



of cabs will overestimate the scattering fraction at high 
energies, and therefore it should not be used when broad 
band energy data is available. Moreover, the column den- 
sity of the warmabs component sets only a lower limit to 
the amount of Compton scattering since there may be 
fully ionised material which cannot be identified via line 
absorption but still contributes to Compton scattering. 



Our final model consisted of disc blackbody 
and blackbody components and one Gaussian 
emission feature at ^6.6 keV, modified by photo- 
electric absorption from neutral and ionised material 
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Fig. 6. EW of the broad emission line with respect to 
the narrow absorption lines of Fexxvi Ka (upper- left), 
Fexxvi K/3 (upper- right), Fexxv Ka (lower- left) and 
Fexxv K/3 (lower-right). 



(tbabs*cabs*warmabs*(diskbb+bbodyrad+gau) , here- 
after Model 2). 

The fits with Model 2 were acceptable, with xl be- 
tween 1.0 and 1.2 for -550-700 d.o.f. for obs 4-9. The 
parameters of the best-fit with this model are given in 
Table U] and the residuals of the fit and unfolded spec- 
tra are shown in Figs. [7] and \E\ For each observation, 
we used the best-fit continuum before adding the broad 
emission line and the photo-ionised absorber as an ionis- 
ing continuum for the latter. The de nsity of the plasm a 
was set to be n = 10 cm~^ following l Ueda et al.l ()2004l) . 
The luminosity was defined in the 0.013-13.6 keV band 
for each spectrum. Interestingly, the value of the ionisa- 
tion parameter for obs 4 in Table |3] does not refiect the 
relative variations in the ionisation of the plasma with 
respect to obs 6-9 that we inferred in Sect. 14.1.11 The 
most likely reason for this is that the model that we are 
using as an ionising continuum for the photoionised ab- 
sorber is not appropriate. If this is true, we expect to find 
the larger discrepancy between the ratio of Fexxv and 
Fexxvi lines and the ionisation parameter in obs 4, since 
this observation shows a column density for the ionised 
plasma larger by more than a factor of two compared to 
the other observations and this will cause severe obscu- 
ration of the continuum. Therefore, we re-fitted again all 
the spectra but using for all the observations the ionising 
continuum of obs 6. We chose obs 6 as the more realistic 
continuum since it has the most highly- ionised plasma and 
consequently more transparent to the incident radiation. 
We obtained log(^) = 4.18 ±0.15, 4.39 ±0.17, 4.33 1°;^^, 
4.00l[^;i^ and 4.09t^;i5 erg cm s^^ for obs 4, 6-9, respec- 
tively. The other parameters and the quality of the fit did 
not change, as expected. These values are consistent with 
the line ratios of different species inferred from Table [3] 
and Fig. 2] and therefore we use these values rather than 
the ones in Table|4]in the remaining sections of this paper. 



Finally, we note that obs 4, 8 and 9 show still after the 
fit small residuals at the energies of '^2.45 and 2.62 keV 
(see Fig. [7]), that we identify with absorption from Sxv 
and S XVI. Since we account for all the absorption fea- 
tures with the self-consistent model warmabs, the exis- 
tence of residuals indicates that the abundance of Sulphur 
is most likely larger than the solar value. Another possi- 
bility would be the presence of a second absorber with a 
lower ionisation degree, but in this case, we should ob- 
serve also residuals for other elements in the plasma. We 
did not attempt to fit the Sulphur abundances due to the 
small significance of the features, but we note that an 
overabundance of the S /Fe ra tio with respect to the solar 
value was also pointed out bv lUeda et al. (2004) based on 
Chandra HETGS data. They argued that this overabun- 
dance could indicate an underestimation of the column 
density of Fe due to the non-inclusion of reemission lines. 
Our data could be showing the same scenario, since we 
also observe small residuals at the energy of the Fexxvi 
edge, which could disappear with an increase of the col- 
umn density of Fe. Alternatively, the remaining residuals 
at the Fexxvi edge in obs 4 and 9 could represent the 
limit of applicability of the warmabs model for high col- 
umn densities (warmabs assumes a uniform ionisation for 
the absorber even if a large column density is present, 
which is not self-consistent). A small absorption feature 
remains also at ~1.88 keV in obs 8 and 9, which could 
be Si XIII, or systematic residuals at the Si edge due to a 
deficient CTI calibration (see Sect. [5]). 

4.1.3. Choice of continuum model 

In previous section we showed that using the same ionis- 
ing continuum for all the observations caused a significant 
change in the value of ^ for obs 4 and 9 and smaller, non- 
significant, changes in the other observations. From this, 
we infer that at least for obs 4 and 9 the continuum may 
not be unique. Continua that result in a similar quality of 
the fit in the 2-10 keV band may yield significantly differ- 
ent ionising continua when they are extrapolated to the 
0.013-13.6 keV energy band. 

There are several reasons for obtaining a "non- 
realistic" continuum in the limited XMM-Newton band- 
pass. First, different combinations of thermal and non- 
thermal components can be used for spectral fitting and 
the quality of the data does not usually allow to discrimi- 
nate among them. Second, if the broad Fe line is caused by 
reprocessed emission, the "reprocessed" continuum should 
be also properly accounted for. A last possibility is that 
changes in the neutral and ionised absorber are mimicking 
a change of continuum, as it is the case in d ipping sources 



( Boirin et al. 2005 : Diaz Trigo et al. 20061) . We examine 



these possibilities in detail below and in next section. 

We substituted the continuum components of the fits 
reported in Table 2] by alternative non-thermal compo- 
nents. We found that the quality of the fits did not 
change. As an example, obs 6 had a xt oi 1.12 for 
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Fig. 7. Obs 4, 6 to 9 are shown from left to right and from up to bottom, (a) 1.7-10 keV EPIC pn (black), 
and 0.8-1.8 keV RGSl (green) and RGS2 (blue) GX13+1 spectra fitted with a disc-blackbody (diskbb), a black- 
body (bbodyrad) and a Gaussian (gau) components modified by absorption from neutral (tbabs) and ionised 
(cabs*warmabs) material, (b) Residuals in units of standard deviation from the above model. 
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Fig. 8. 0.8-10 keV unfolded spectra. Solid, dashed, dot-dashed and dot-dot-dashed lines represent the contribution of 
the total model, disc, blackbody and Gaussian components, respectively. 



701 d.o.f. when fitted with Model 2. Substituting the 
blackbody component by a power-law component in 
Model 2 (tbabs*cabs*warmabs*(diskbb+po+gau) , here- 
after Model 3) we obtained xl = 1-12 for 701 d.o.f.. 
The photon index, F, and normalisation of the power- 



law were 2.2 ± 0.3 and 1.1 ± 0.4, respectively, and the tem- 
perature and normalisation of the disc blackbody compo- 
nent 1.85 to 06 ^^'^ 461^4, respectively. Substituting 
the disc blackbody component by a power-law component 
in Model 2 (tbabs*cabs*warmabs* (bbodyrad+po+gau) , 
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Table 4. Best-fits to the 0.8-1.8 keV RGS and 1.7-10 keV EPIC pn persistent spectra for all the observations using 
the tbabs*cabs*warmabs* (diskbb+bbrad+gau) model. The column density of the cabs component was tied to the 
column density of the photo- ionised absorber (see text), khh, ^dbb and fcgau are the normalizations of the blackbody 
component, disc blackbody and Gaussian emission feature, respectively. fcTbb and kT^hh are the temperatures of the 
blackbody and disc blackbody, respectively. Eg^u and a represent the energy and width of the Gaussian feature. TVfj'^ 
and TVy'^''™'^'"' are the column densities for the neutral and ionised absorbers, respectively. a^, and v are the ionisation 
parameter (in units of erg cm s~^), the turbulent velocity broadening, and the average systematic velocity shift of the 
absorber (negative values indicate blueshifts). Ftot is the unabsorbed total flux, and Fgau is the unabsorbed emission 
Gaussian flux. The width (cr) of the Gaussian emission line gau was constrained to be <1 keV in the fits. 



Observation No. 



Parameter 

fcTdbb (keV) 

fcdbb [(Rin/Dio)^ cos6l] 

fcTbb (keV) 

fcbb [(R»n/Dlo)'] 

Sgau (keV) 
a (keV) 

EW (eV) 



Comp. 
dbb 

bb 

gaui 

) 

tbabs 



Nl^'' (lO^^ cm-2) 

log(C) (erg cm s"^) 
(Tv (km s~^) 
V (km s"""") 



warmabs 



0.86 

492 t«« 



1.24 ±0.02 



310 



+37 



6.56 

0.88 
0.023 ±0.006 
299 ± 77 



2.86 ±0.05 



r ' _g.g 
4 41 +0-11 

^•^^ -0.19 

-165 
-75 
-H420 
-360 



580 : 
-2310 



1 Ifi -1-0.19 
^•J^O -0.07 

226 til 



1 cc-l-0.24 
l-66_o.03 



74 



+8 
-35 



6.71 



-1-0.12 
L -0.16 
f^ 77+0.23? 
'J-' '-0.15 

0.0081°°°^ 



6-6 ta.? 

4.39 ±0.17 

33Q0 +l''85 
ooau _22io 

-3690 



-0.08 

2011:11 



1 fiO+°°* 

r.DU _o.02 



81 



+7 
-31 



6.63 ±0.11 

n 79 +0.08 
U. -0.12 

U.UU< _o.ooi 

88 tit 



6.6: 

4.34: 
2130: 
-2160 



-2.9 

-0.4 

-0.22 

-0.08 

-1220 

-920 

+270 

-180 



i.8olg:li 

38±8 



-I 91 +0.08 
J^-^J^-0.10 

3+34 
-38 



113: 



6.59 

0.78 tVol 
0.014 

19511° 



2.86 ±0.02 2.92 ±0.02 2.77 ±0.05 



6.0: 



-1+6.6 
-2.0 
3^97+0.19 



0.06 

+490 
175 



480: 
-2130 ±480 



„ 77+0.06 
U. ( ( -0.04 
f 188 
-110 



622: 



1.28 ±0.01 
272 til 

6.69 ±0.06 

n 77+0.06 
' ' -0.08 

0.017 ±0.003 
235 ± 42 



9 07 +0.06 
^.»' -0.04 



7 -r +0.5 
■ — 2.3 

4.23 ±0.10 
975 

^'^ -710 

-2400 ± 360 



x'(d.o.f.) 


1.11 (586) 


1.12 (701) 


1.04 (686) 


1.09 (548) 


1.23 (575) 




F^Jio keV (10-°' erg cm'^ s"!) 


9.6 


9.5 


8.7 


8.0 


8.3 


^^1-10 keV (10-" erg cm-2 s'l) 


24.3 


7.9 


7.1 


14.7 


17.6 


Exposure (ks) 


5.2 


8.9 


13.5 


6.5 


6.7 



hereafter Model 4) we obtained xt = 1.14 for 701 d.o.f.. 
The photon index, F, and normalisation of the power-law 
were 2.27 ±0.04 and 3.1 ±0.2, respectively, and the tem- 
perature and normalisation of the blackbody component 
1.12 ±0.02 keV and 28119, respectively. We did not find 
significant changes in the parameters of the narrow fea- 
tures when using different continua. In contrast, the EW 
of the Gaussian component increased from 91 123 in 
Model 2 to 201 ± 41 and 348 tfe eV in Models 3 and 4, re- 
spectively. Finally, substituting the power-law component 
in Models 3 or 4 by a cutoff or broken power-law compo- 
nents gave an unconstrained cutoff or break and a similar 
fit to the power-law. 

We also obtained fits of the same quality when using 
different abundances. However, we obtained a change in 
the EW and breadth of the Gaussian component. As an 
example, in obs 6 the breadth and EW of the broad emis- 



sion line increased from 91_23 lo 129 _2i eV when we 



substituted the s olar abundances of Anders fc Grevess3 
(fl989) by those of IWilms et al] (|2000l) in Model 2. 

In summary, we found that the parameters of the warm 
absorber were stable to different choices of continuum 
components and abundances. In contrast, the parameters 
of the broad Fe line changed significantly, the main rea- 
son being that the quality of the fit before introducing 
the broad Fe line was worse for Models 3 and 4 than for 
Model 2. 



4.1.4. Modeling of the broad Fe line and continuum 
variability 

We examined the possibility that spectral changes be- 
tween different observations are due to variations of the 
warm absorber and the reprocessed component instead of 
variations of the continuum. A motivation for this is Fig.[H 
The resemblance of this figure with Fig. 4 of iBoirin et al 
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( 2005[) indicates that obs 8-9 (7) could be deep (shallow) 
dipping states with respect to obs 6 (see Sect. I3.ip . If this 
is true we expect the continuum to be relatively stable and 
the dipping states to have a less ionised absorber and with 
a larger column density compared to the persistent ones 
(jBoirin et alj 120051: iDfaz Trigo et al.ll2006l ). Indeed, obs 8 
and 9 show a less ionised absorber with visible features of 
Fexxv, compared to obs 6 and 7. Similarly, obs 1-3 show 
a less ionised plasma than obs 4. 

In order to test the continuum variability, we 
first substituted t he Gaussian c ompo nent by the 
reflionx model ( Ross k. FabianI l2005l) to account 
for reflection by an ionised, optically thick, illumi- 
nated atmosphere of constant density. Since reflionx 
only provides the reflected emission, we added also a 
power-law component to account for the direct emission 
(tbabs*cabs*warmabs* (diskbb+bbodyrad+po+ref lionx) , 
hereafter Model 5), and we coupled the photon index of 
the illuminating power-law component to the index of 
reflionx. If the line is produced by reprocessed emission. 
Model 5 will yield a more realistic continuum, necessary 
to perform a simultaneous fit to several observations. 
We note that while reflionx is the most adequate of 
the available public models for our purposes since it 
includes self-consistently line emission and accounts for 
ionised reflection, n iore recent refined mode ls such as 
those presented in Garcia fc Kallman ( 2010l ) point to 
important differences in the resultant spectrum when e.g. 
new atomic data are considered (see Sect. 15. 2| ). 

For the simultaneous fit we considered obs 6-9 only. 
We left out obs 4 since it is in a different track in Fig. [2j 
together with obs 1-3 and this may indicate already sig- 
nificant changes in the continuum (a significant increase 
in Compton scattering due to a larger column density of 
plasma could be also the reason for the different track of 
obs 1-4). We imposed the same continuum for obs 6-9 
(including index and normalisation of the power-law com- 
ponent) and allowed only changes in the neutral absorber 
and the ionisation degree and normalisation of the reflec- 
tion component. We kept the parameters of the warm ab- 
sorber fixed to the values shown in Table H] but with the 
values of the ionisation parameter obtained when the ion- 
ising continuum of obs 6 was used to fit all observations 
(see Sect. 14.1^ . Finally, we allowed a global difference in 
normalisation among observations to account for a poten- 
tial different amount of Compton scattering. 

We obtained a good quality for the simultaneous fit of 
obs 6-9 {xl = 1.21 for 2537 d.o.f.). Table[5]shows the best- 
fit parameters and the residuals of the fit and unfolded 
spectra are shown in Fig. [Sj 

We could model the emission line with the reflection 
component and, as expected, the reprocessed component 
contributes significantly to the continuum emission. Thus, 
the interpretation of the continuum parameters should be 
done with caution. Interestingly, the large difference in the 
curvature of the spectrum at 7-10 keV between obs 6-7 
and 8-9 can be well reproduced by a difference in the 
ionisation and normalisation of the reflected component. 



In contrast, we were not able to reproduce the broad iron 
line with reflection models which consider a blackbody as 
an incident spectrum since they predict EWs that are too 
small compared to the observed lines. 

Despite the success of Model 5 to reproduce the various 
spectral features of obs 6-9 while assuming a constant 
continuum, there are still some systematic residuals in the 
best- fit. The residuals at the absorption edge at the energy 
of Fe XXVI in obs 9 were already present in individual fits 
(see Sect. I4.1.2[) . However, the simultaneous fits show a 
slight curvature in the residuals which is not present in 
the individual fits and could be due to shortcomings of the 
reflection component, a different origin for the emission 
line, or an indication for continuum variability. Finally, 
if we attribute the differences in the global normalisation 
constant to Compton scattering of fully ionised plasma, 
the ~30% difference in normalisation between obs 6-7 and 
8-9 implies that obs 6-7 are observed through an extra 
column density of 10^'^ cm~^ of fully ionised plasma. 

4.2. RGS spectral analysis 

We examined the 0.8-1.8 keV (9.9-22.3 A) first order RGS 
spectra to constrain the A^h in the direction of the source 
and to search for the signature of narrow absorption and 
emission features. We note that the second order spectra 
have too few counts to add any value to the fits. 

We could fit the RGS spectra of all the observations 
with a continuum consisting of a disc blackbody modi- 
fied by photo-electric absorption from neutral material. 
We obtained a C-statistic of between 450 and 571 for 
-434-436 d.o.f. for obs 4 to 9). The values of A^h in 
units of 10^2 cm-2 were 2.72 ± 0.06, 2.73 ± 0.05, 2.80 
2.80 ±0.06, 2.76 ±0.07 and 2.75 tJJ;!^^ for obs 4 to 9, re- 
spectively. None of the observations showed significant 
narrow features. 

To evaluate the effect of the chosen continuum in 
the value of iVn, we substituted the disc blackbody by 
a power law. We obtained a C-statistic of between 448 
and 570 for ~434-436 d.o.f. for obs 4 to 9. The values 
of TVh in units of lO^^ cm'^ were 2.53 ± 0.25, 2.53 ± 0.19, 
2.80 ± 0.19, 2.71 ±0.20, 2.41 ±0.26 and 2.16 1^;^^ for obs 4 
to 9, respectively. 

Finally, to evaluate the effect of the absorption com- 
ponent in our fit, we substituted the tbabs component by 
the more recent tbnew component and re-fittcd the com- 
bined spectra. The C-statistic and values of A^h did not 
change. 

In summary, we obtained consistent values of A^h for 
different models of interstellar absorption but significantly 
different values of A^h when using different continuum 
components. This is expected due to the combination of 
limited energy band of the RGS and low statistics of these 
observations. Therefore, we conclude that the RGS spectra 
alone cannot constrain the value of A^h for these observa- 
tions due to the poor statistics and the limited bandwidth, 
necessary to determine the continuum. 
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Fig. 9. Left: (a) 1.7-10 keV EPIC pn, and 0.8-1.8 keV RGSl and RGS2 GX13+1 spectra fitted with Model 5 (see 
text), (b) Residuals in units of standard deviation from the above model. Right: 1-10 keV unfolded spectra. Solid, 
dot-dot-dashed and dot-dashed and dashed lines represent the contribution of the total model, disc, blackbody and 
reflection components, respectively. Obs 6-9 are shown in black, blue, green and red, respectively. 

Table 5. Best-fits to the 0.8-1.8 keV RGS and 1.7-10 keV EPIC pn persistent spectra for all the observations using 
the const*tbabs*cabs*warmabs* (diskbb+bbrad+po+ref lionx) model (see caption of Table 0] for definitions). T is 
the index of the incident power-law and is coupled to the index of ref lionx. kpo is the normalisation of the power-law 
component and ^ and krefUonx are the ionisation parameter and the normalisation of the reflection component. The 
parameters of warmabs and cabs have been fixed to the values of Tabled 

Observation No. 6 7 8 9 



Parameter 
fcTdbb (keV) 

fcdbb [(Rin/Dlo)2 cose] 

fcTbb (keV) 

fcbb [(R»n/Dlo)^] 



kpo (ph keV"^cm"^s"^ at 1 keV) 
log(f) (erg cm s"^) 

^ref lionx (10 ) 

N'^" (10^2 cm-2) 
Const 



Comp. 
dbb 

bb 

po 

ref lionx 

tbabs 

const 



3.94 in? 

o 1 +0.6 
■^•^ -0.3 



0.86 IH^e 

I +27 



189: 



1.23 
199 



+0.01 
-0.04 

[-25 



1.94 ±0.02 
< 0.006 



q oq +0.05 
O.OO _o.08 
q 9 +0.5 
■^•^ -0.3 



3.41 ±0.04 

q o +0.5 
■^-O -0.3 



o on +0.01 
O.OU _o.o4 

4 6+0-4 

-0.2 



3.13±0.03 3.17±0.03 3.17±0.03 2.94±0.03 



1 (f) 



1.02 ±0.03 1.29 ±0.03 1.37 



x'(d.o.f.; 



1.21 (2537) 



In order to increase the sensitivity to narrow lines and 
edges, we combined the RGS spectra of several observa- 
tions. We added the RGSl and RGS2 spectra of obs 4, 
8 and 9 (set 1) and separately of obs 5, 6 and 7 (set 
2). The reason to distinguish between sets is that the 
former set of observations shows possible "obscuration" 
events (see Sect.|4]) and, if confirmed, such events indicate 



a plasma with lower ionisation degree than the second set 
of observations. We could fit the resultant spectra with 
the same model as the individual observations. We ob- 
tained a C-statistic(d.o.f.) of 318(280) and 390(280) and 
iVH=2.74±0.04 and 2.78j:!^:^4 xlO^^cm'^ for sets 1 and 
2, respectively, when using a disc blackbody continuum 
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and of 292(280) and 385(280) and A^h=2.33 ± 0.15 and 
2.67±0.11 xlO^^cm^^ when using a power law. 

None of the merged data sets shows any additional 
narrow features and including absorption from a photo- 
ionised plasma in the model does not improve the fit qual- 
ity. We calculated upper limits for the existence of the 
Mg XII absorption feature de tected in the Cha ndra obser- 
vation with an EW of 2 eV ( Ueda et al]l2004l ). Assuming 
that the feature is narrow (width fixed to 0), we obtained 
upper limits to the EW between 1 and 3 eV for obs 4-8 
for a feature with no velocity shift. For obs 9 we obtain 
a non-significant detection of an absorption feature with 
EW = 2.5tl;l eV. Then, we assumed velocity blueshifts 
up to 4000 km/s in steps of 500 km/s and calculated up- 
per limits at the corresponding wavelengths. We obtained 
upper limits to the EW between 0.5 and 4.2 eV for obs 4, 
7-9. For obs 6 we obtain a non-significant detection of an 
absorption feature at a blueshift of 2000 km/s with an 
EW = 2.0tli eV. These results are consistent with the 
predicted EW for Mgxii, <;i eV, for an ionised absorber 
of log 4, as fitted in the EPIC pn exposures. 

4.3. OM data 

We extracted images and light curves from the U, UVWl 
and UVM2 filter OM exposures for obs 4-9. 

The images show a source in the U filter exposures con- 
sistent with the position of GX 13-1-1, which is strongest in 
obs 4 and weakest in obs 7. The source was not detected in 
the UVWl or UVM2 filters, except for obs 5 where there 
is a detection with the UVWl filter with a magnitude 
of 20.7 ±0.7. For obs 4, 5, 7 and 8 we obtained average 
U optical magnitudes of 19.8 ±0.4, 19.7 ±0.2, 20.8 ±0.7 
and 20.5 ± 0.5, respectively. We note that GX 13+1 is not 
detected in the U exposures in obs 6 and 9. 

Obs 4 showed a significant variation in the U magni- 
tude with time, from 20.4 ±0.4 to 19.2±0.1 within a 9 ks 
interval. Unfortunately, we cannot determine the presence 
of a modulation with orbital phase due to the few detec- 
tions of the source. Therefore we do not use the OM data 
any further. 



fit. However, we found a significant change in the param- 
eters of the broad emission line when such non-thermal 
components were used (see Sect. 14.11 and l5.2|) . 



Lin et al.l (|2007l ) showed that the commonly used mod- 



els for thermal emission plus Comptonisation were not 
favoured for the "Soft State" of two atoll sources because 
they failed to satisfy various desirability criteria, such as 
the Lx oc T^ evolution for the multi-colour disc blackbody 
component and similarity to black holes for correlated tim- 
ing/spectral behaviour. Similarly, in a systemati c analysis 



of 16 NS LMXBs observed with XMM-Newton, iNg et al. 
( 2010l ) found that a continuum consisting of a disc black- 
body and blackbody components was favoured in 80% of 
the cases. In this model, the blackbody component was in- 
terpreted as emission from a radiation pressure supported 
boundary layer. The data analysed in this work give fur- 
ther support to a continuum with two thermal compo- 
nents. However, we note that since XMM-Newton effec- 
tive area is already very low at ^10 keV, we cannot de- 
cide based on this data whether an additional power-law 
component extending to high energie s, like the one de - 



tected by INTEGRAL up to ~40 keV (jPaizis et al.ll2006D 



would be requested above 10 keV. If confirmed, the pres- 
ence of a reprocessed component representing reflection 
by an ionised, optically thick, atmosphere illuminated by 
a power-law component (see Sect. 14.1"!^ gives indirect ev- 
idence of the existence of a power-law component extend- 
ing to high energies. 



Ueda et al.l (|2004l ) used a blackbody and disc black- 
body to fit the 3-25 keV RXTE/PCA spectrum of 
GX 13+1 and obtained a satisfactory fit but the values for 
the temperatures of the disc and blackbody components 
were significantly different than those reported in Table SI 
Taking into account the flux and the ionisation of the disc 
wind, the RXTE/PCA obse rvation shou l d be in a simi- 
lar state to obs 7. However, Ueda et al. ( 20041) report a 



temperature of 1.52 ± 0.03 keV and of 2.56 ± 0.02 keV for 
the disc and blackbody components, much larger than the 
ones found in this work for obs 7. Clearly, the absolute pa- 
rameters of the continuum should be taken with caution 
in fits with limited energy band coverage. 



5. Discussion 



5.1. Photoionised absorption 



We analysed five XMM-Newton observations of the 
LMXB GX 13+1 taken during 2008 to investigate the vari- 
ability and origin of the disc wind present in this source. 

The 0.6-10 keV EPIC pn lightcurves were highly vari- 
able with two observations showing energy- dependent ob- 
scuration, similar to the classical dipping sources. 

The X-ray continua of the "least variable" intervals 
were well fitted by a model consisting of a blackbody 
and disc-blackbody components absorbed by neutral and 
ionised material. We could substitute the blackbody or 
disc-blackbody components by a non-thermal component 
such as a (cutoff) power-law or a Comptonisation com- 
ponent without significantly worsening the quality of the 



Absorption from highly-ionised species was present in all 
the observations of GX13+1 analysed in this work. 

An absorption featu re consisten t with Fc xxvi was first 



discovered by ASCA (jUeda et al.M2001.) . XMM-Newton 



observations in 2000 (obs 1-3) showed a complex of ab- 
sorption lines and edges consi stent with absorption from 
Caxx, Fexxv and Fexxvi (Si doh et al.li2002,'). similar to 



the complex found in obs 4, 8 and 9. lUeda et al.l (|2004f ) 



found additional features of Sxvi, Sixiv, Mgxii, Mnxxv, 
Cr XXIV and Ar xvii in a posterior Chandra HETGS ob- 
servation in 2004 and measured a blueshift of ^^400 km 
and an ionisation parameter of log(^) ^ 4.1-4.7 for the 
plasma responsible for the absorption. A plasma density 



Di'az Trigo et al.: GX13+1 



15 



^10 cm"'^ and a launching radius for the wind of ^10 - 
10^^ cm was inferred. 

The depth of the absorption features in the XMM- 
Newton 2008 observations changed significantly in 
timescales of few days, and m ore subtle va,riatio ns are seen 
in shorter, hours, timescales ( Sidoli et al. |2003) . 

We found column densities for the absorber between 6 
and 17 xlO^^ cm~^ and ionisation parameters, log(^), be- 
tween 4 and 4.4 erg cm s~^. Obs 4 shows a column density 
larger by more than a factor of two with respect to the 
other observations, and which is poorly constrained. This, 
and the well constrained velocity indicate that the lines 
lie in the saturated region of the curve of growth, where 
the EWs of the lines do not depend on the column den- 
sity anymore but increase with the velocity broadening. 
We observed a correlation between the ionisation param- 
eter and the turbulent velocity of the warm absorber (see 
Fig.Hni). 
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Fig. 10. Ionisation parameter of the warm absorber with 
respect to its turbulent velocity. 

We found blueshifts between --2100 and 3700 km s"^ 
for obs 4-9. Such blueshifts seem too large, compared to 
the shifts of -400 km s^^ obtained bv lUeda etaP (|2004l) 
for GX 13-1-1 and of —300-1600 km s~^ in microquasars 
'e.g. Kallman et al. 2009t Ueda et all l2009i see however 



Miller et al.l ( 201 ll ) for a recently reported blueshift of 
3100 ±400 km/s in a Chandra HETGS observation of 
IGR J17480-2446). Unfortunately, we do not detect any 
discrete absorption features from the warm absorber in 
the RGS, which has a much better energy resolution than 
EPIC. The lack of spectral features in the RGS wavelength 
range is consistent with the high degree of ionisation of 
the absorber (see Sect. 14. 2p . Thus, a more accurate mea- 
surement of the blueshifts is not possible for these obser- 
vations. However, it seems plausible that the shifts de- 
tected with the pn camera are overestimated compared to 
the real blueshifts due to residual calibration uncertainties 
(see Sect. liXTj) . 

We can calculate the distance between the ionising 
source and the slab, r. Since ^ = L/uer^ and Ue — 



^warmabs^^ (wlicrc d represents the thickness of the slab 
of ionised absorbing material), we can calculate r as 
{L/^ 7Vg^™'*'=")(d/r). Considering plausible values of d/r 
range between 0.1 and 1, we obtain values of r between 
3x10^ cm and 1x10^^ cm. Further, assuming that the 
value of r should not change significantly between ob- 
servations, we estimate r as 1-3 x 10 ^° cm, in excellent 



agreement with the value obtained by lUeda et al.l ()2004l ) 



of 2 X 10^° cm. 

The characteristics of the wind allow us to discrimi- 
nate among the possible launching mechanisms, namely 
thermal, radiative, magnetic or a combination of them. 
Thermal driving is effective at large distances from the 
central compact object, black hole or NS, where the ther- 
mal velocity exceeds the local escape velocity. Thermally 
driven models are thus most effective in producing slow 
winds and at large radii. The distance at which the wind 
is launched in GX 13-1-1, 1 -3 x 10^° cm , point s to a ther- 
mal origin, as discussed bv lUeda et al. (2004). Moreover, 
the fact that t he only known NS LMXBs showing winds 
are GX13-H1 (lUeda et al.l l2004h and IGR J17480-2446 
(jMiller et al.ll201ll ) is a logical consequence of the wind 
launching mechanism being thermal. As X-rays heat 
low density gas to a temperature T~10'' K, the sur- 
face of a heated disc is expected to either puff up and 
form a static corona or emit a thermal wind, depend- 
ing whether t he thermal velocity exceeds the local es- 
cape velocity (Begelman et al.l 1983t Woods et al. 1996 



Proga fc KallmanI I2OO2I) . IWoods et al.l (|l996[ ) estimated 
that for a radius of ~2.6 x 10^° cm, an isothermal wind 



would develop for luminosities above —3x10'^^ erg s" 
while such a wind would be inhibited by gravity at lower 
luminosities (see their Fig. 17). For larger radii, a steadily 
heated free wind coul d develop already at lumi nosities be- 



low-3xl03'^ ergs~^ |Proga fc KallmanI (|2002[) found that 



including the radiation force due to electron scattering, 
the radius at which the maximum of mass flux density oc- 
curred could be reduced by a factor of 20 compared to the 
calculations of IWoods et al.l (|l996l ). Taking these values 
into account, the luminosity, and not the radius, seems 
to be the critical parameter for the existence of a wind 
in the classical dipping sources. Indeed, all the dipping 
sources showing ionised absorption for which no signifi- 
cant blueshifts have been detected (Cottamet al.ll2001 



'Sidoh et al."200l':'juett fc Chakrabartvl2006HBoirin et al 
,2005i : ilaria et al. 2007) have a luminosity below the one 
needed for the wii id to overcome gravity in these mod- 
els (see Table 1 of Di'az Tr igo et al. '2006) , with the pos- 
sible exception of MXB 1658—298, which has a luminos- 
ity close to the critical value. In contrast, X 1624—490, 
GX13-hl and IGR J17480-2446 have luminosities of 
4.7 X 10^'', 6-9 X 10^"^ and 3.7 x 10^'^ erg/s, respectively, 
close to or above the critical luminosity. From these 
sources, two show clearly outflowing winds (see above) and 
the third one, X 1624—490, shows a significant blueshift of 
607 

^342 km/s for the most highly ionised ab sorbing com 



pone nt in Chandra HETGS observations (IXiang et al 
l2009i) . In summary, a thermal wind launching mechanism 
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can explain the appearance of hot static atmospheres or 
slow disc winds in all high inclination NS LMXBs and its 
absence in low inclination NS LMXBs. We caution that 
Compton scattering in a hot atmosphere or disc wind with 
column densities above 10^^ cm~^ can reduce significantly 
the flux of the central source seen by a distant observer. 
Therefore, a potential larger luminosity than "observed" 
(as the luminosity of the system could be significantly un- 
derestimated due to Compton scattering) and the radi- 
ation force due to electron scattering when a large col- 
umn density of highly ionised plasma is present, should 
be always considered before ruling out the thermal wind 
launching mechanism. 

Finally, if the thermal wind launching mechanism is in 
place for all NS LMXBs, we expect that a disc wind will 
develop for systems brighter than GX 13-1-1, known as Z- 
sources , si nce their luminos ity is well above the critical 
luminosity. lUeda et al. estimated the mass outflow 

rate in GX13-hl to be ^ 0.7-2.7 x lO^^ g s'^. They as- 
sumed that the plasma had an equatorial geometry and 
was outflowing from both sides of the disc to define a solid 
angle of ^I/Att = 0A. In our observations, taking a similar 
solid angle into account we calculate outflow rates of 0.1- 
0.2 X 10^^ g s~^ for an outflo w velocity of 400 km s^^. 



as found bv lUeda et all (|2004[) . and of 0.5-1 x 10^" g s"^ 
for an outflow velocity of 2000-3500 km s^^, as found in 
these observations. This outflow rate is comparable with 
the mass accretion rate of 10"'^^ g s"-^, as already pointed 



out by lUeda et"aL ( 2004 ) , and implies that the disc wind 
could play an important role in the wh ole dynamics of the 
accret ion disc for such bright sources. Proea fc Kallmaiil 
(I2OO2I) calculated the mass-outflow rate for a NS accreting 
at a rate of 1.3 x 10^^ g/s and found mass-outflow rates be- 
tween 3.9 x 10^^ and 0.65 x 10^® g/s for different values of 
X-ray attenuation (see their Table 1). Unfortunately, from 
their simulations we cannot estimate what is the effect of 
increasing the mass accretion rate in the mass-outflow rate 
for a given value of X-ray attenuation. However, the fact 
that for the same attenuation, an increase in accretion 
rate results in an increase of column density of the ionised 
plasma (see e.g. models L2 and H2 in their Table 1) may 
indicate that the mass-loss rate will be larger for sources 
with higher accretion rate, since the radiation pressure due 
to electron scattering will also increase for such sources. 

Although narrow absorption features are not detected 
in the Z-sources, most likely due to a lower inclination or 
to fully ionised winds, if we consider other signatures of a 
disc wind such as broad iron line emission (see Sect. 15.31) 
a study of the strength of the wind in such sources may 
become possible and of utmost importance to evaluate the 
balance between inflowing and outflowing matter in such 
systems. 

5.2. Broad iron line emission 

We found a broad, a ~ 0.7-0.9 keV, emission line in all the 
observations analysed in this work. The energy centroid 



varied between 6.55 and 6.7 keV, consistent with emission 
of highly ionised species of iron (Fe XXi-Fe xxv) . The EW 
of the line increased from ^ 100 eV in obs 6 and 7 to 
-300 eV in obs 4. 

The absolute breadth and EW of the lines was strongly 
dependent on the continuum model. As an example, the 
EW of the line in obs 6 increased from 91 ^23 to 201 ± 41 
(348 i^g) eV when we substituted the disc blackbody 
(blackbody) component by a power-law component in 
Model 2 (see Sect. 14. iT^ . For this highly absorbed source 
we also found a dependency of the breadth and EW of 
the line on the abundances used, which caused a significant 
change in the continuum parameters, especially in the soft 
band. For obs 6, we ob tained an EW of 9 1 1,^ e V using 
solar abundances from Anders &: Grevesse (1 19891) an d of 
129^21 eV with abundances from lWilms et all (|2000[ ). In 
general, continuum models that yielded a higher xt for the 
fit before adding the line resulted in broader lines, since 
the Gaussian component absorbed the deficiencies of the 
continuum model. Hence, it is challenging to investigate 
the origin of the broad iron line based on its absolute pa- 
rameters. In contrast, we can examine the relative changes 
of the line among observations, since any systematic effect 
caused by our choice of continuum, abundances or spectral 
analysis should be the same for all observations. Further, 
since we are comparing observations of one source and 
taken with the same instrument, residual calibration ef- 
fects or properties of the source such as inclination should 
not affect our conclusions. 

The presence of broa d Fe lines in NS LMXBs has been 



extensively studied (see IWhite et al.l 119861 : iHirano et al 
1987tlAsai et al.ll2000l: INe et allbOloT for systematic stud- 



ies with different X-ray observatories). In these systems, 
broad lines could arise in the inner accretion disc by flu- 
orescence following illumination by an external source of 
X-rays (e.g. the boundary layer where the accretion disc 
meets the star), and be broad ened by relativistic effects 
near the compact object (e.g. Revnolds fc Nowa"^ 12003 : 
Fabian fc Miniuttil [200I iMattI 120061 ) . Ahernativelv. they 
could originate in the inner part of the so-called accre- 
tion disc corona, formed by evaporation of the outer lay- 
ers of the disc il luminated by the emission of the cen- 
tral object (e.g., IWhite fc Hdtlll982l: iKallman fc Whitel 
1989t ) and be broadened predominantly by Compton scat- 
terin g fe.g. IPozdnvakov et al.lll979l:ISunvaev fc Titarchuk 
1980l ). A third possibility, is that they originate in a 
partially ionised wind as a result of illumination by 
the central source continuum photons and are broad- 
ened by electron downscatte ring in the wind environment 
(jLaurent fc TitarchukI 120071 ). The first interpretation dif- 
fers from the other two in that Compton scattering does 
not suffice to broaden the lines to the observed values, 
hence relativistic broadening has to be invoked. Therefore, 
to discriminate between the different interpretations it is 
key to determine the characteristics of the plasma where 
the lines originate, mainly its temperature and state of 
ionisation, by comparison with detailed reflection models. 
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Recently, Garcia fc Kallman ( 20101 ) presented new 
models for illuminated accretion discs and their structure 
implementing state-of-the-art atomic data for the isonu- 
clear sequences of iron and oxygen and increasing the en- 
ergy, spatial and angular resolution compared to previous 
works. They found EWs of the Fe K emission line be- 
tween 400 and 800 eV for the cases with log(^) between 
1 and 3. The line EWs decreased rapidly for higher val- 
ues of with EW ~ 40 eV for log(0 = 3.8. The EWs 
of the lines found in this work range between ~90 and 
~300 eV, indicating a plas ma with log(£) bet ween ^^3.7 
and ^ (see also Fig. 2 of Garcia et al. 2011 ). in broad 
agreement with the values of log(^) ~ 3.3-3.9 found for 
the reflection component used to model the emission line 
in Sect. [ 4.1.41 For t hese values of the ionisation parameter 
Garcia et al.l ( 2011 ) found that the most prominent lines 
had centroids between 6.65 and 6.74 keV, slightly larger 
than the values found in this work, 6.56-6.71 keV, but still 
consistent within the errors. 

In summary, we found that the lines observed in this 
work can be explained by ionised reflection of the disc. In 
this model, the breadth of the lines results primarily from 
the Compton scattering in the hot atmosphere or corona, 
and relativistic effects are not required as a broadening 
mechanism. 

However, the contemporaneous evolution of warm ab- 
sorption and broad Fe line emission for the observations 
analysed in this work indicates that they probably form 
in the same environment, n amely the disc w ind (see 
sect. 15. 3p . As pointed out by ISim et all ()2010al ). in the 
presence of a disc wind we expect a scattered/reprocessed 
component that is qualitatively similar to that produced 
by standard disc reflection models. In particular, a moder- 
ately strong Fe line is present in both models. In contrast, 
it is difficult to justify the contemporaneous changes of the 
broad iron line and the absorption features if the former 
is produced at the inner disc or the boundary layer, as 
expected in the classical reflection scenarios, and the lat- 
ter at the outer disc, as shown in Sect. 15.11 We note that 
strong broad Fe lines were predicted for LMXBs which 
had a wind with a Thomson scattering optical depth close 
to unity far from the compac t object and which is illum i- 
nated by a hard component (jLaming &: TitarchukI 120041 ). 
In that model, the optical depth in the Fe K continuum 
is about 1-3 times higher than that due to electron scat- 
tering. The wind found in obs 4, for which the Fe line is 
strongest, has a column density of 2 x 10^"^ cm^^, but this 
may be only a lower limit to the real column density. The 
fact that the absorption lines are saturated and that the 
absorption edges are not properly fitted points to a failure 
of the warmabs model when the optically thick limit is ap- 
proached. Therefore the conditions necessary for a wind 
origin of the broad line are most likely fulfilled. 



Finally, while we cannot assure that a hard 
illuminating compo n ent, as the one assumed by 
Laming fc TitarchukI (j2004l ). is present in this data due 



^ We note that these values of log(^) cannot be directly com- 
p ared with those of th e warm absorber, since the models shown 
in lGarcfa fc KallmanI l|2010l ) have an ionising continuum con- 
sistent of a power law with index 2. 



to the limited energy coverage of XMM-Newton, if such 
component is absent, disc reflection would not be an al- 
ternative mechanism, since strong lines cannot be pro- 
duced with a so ft illuminating component (see Fig. 3 of 



Banantvnell2004[ ). 



Therefore, we favour reprocessing in a hot atmosphere 
or disc wind as the origin for the broad Fe line in NS 
LMXBs and justify this in detail in next section. 

5.3. A disc wind as a self-consistent scenario for 
photoionised absorption and broad iron line 
emission 



iNeil sen fc Lei reported the presence of a broad 

emission line in the faint, hard states and narrow absorp- 
tion lines in the bright, soft states of the microquasar 
GRS 1915-1-105. Since the hard states exhibit prominent 
radio jets, which are absent during the soft states, they 
argued that the broad emission line arises when the jet il- 
luminates the inner accretion disc and the absorption lines 
originate when the powerful radiation field around the 
black hole drives a hot wind off the accretion disc. Since 
GX 13+1 sho ws also jets, which a re associated with hard 
states ( Grind lav fc Seaauis il986l) . we could expect a sim- 
ilar behaviour regarding the appearance of broad emission 
lines and narrow absorption lines as for GRS 1915-1-105. 

In contrast, we observe contemporaneously broad 
emission line and narrow absorption lines similarly to most 
of the classical dippers, and the emission line becomes 
stronger as the hard flux decreases. Moreover, we find that 
the EW and the width of the broad Fe line increase with 
the column density of the absorber (see Fig. 1111 upper pan- 
els) , a decrease of the EW of the broad line with the in- 
crease of the ionisation parameter of the absorber fFig.fTTl 
lower-left panel) and a simultaneous increase in the ionisa- 
tion of the absorber and the energy centroid of the broad 
line (Fig. 1111 lower-right panel) . We can naturally explain 
all these correlations in the frame of a disc wind scenario 
in which both the absorption and the emission are pro 



duced in the wind (Laurent fc Titarchuk 2007: Sim et al 



2010a '). In this scenario, we observe absorption features as 



we look through the wind (predominantly in high inclina- 
tion sources) , while the radiation scattered in the outflow 
produces the broad iron line emission and other distinct 
features such as a Compton hump, which will be more or 
less visible depending on their contribution with respect 
to the incident radiation. Higher column densities make 
scattered/reprocessed radiation more dominant and so we 
observe a most prominent iron line whenever the column 
density of the absorbin g plasma is larger and the ionisa- 
tion smaller. We refer to Sim et al. I (l2010al lbl) for a detailed 
explanation of the total spectra expected at different in- 
clinations and for different characteristics of the wind and 
summarise the relevant results in what follows. 
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Fig. 11. Upper panels: EW (upper- left) and width (upper-right) of the broad iron hue with respect to the column 
density of the warm absorber. Lower panels: EW (lower-left) and energy (lower-right) of the broad iron line with 
respect to the ionisation of the warm absorber. 



At low, ^45°, inclinations the wind does not obscure 
the X-ray source. However, a significant component of 
scattered/reprocessed radiation in the outflow is present 
in addition to the direct emission and responsible for the 
broad iron line emission. The emission is predominantly 
formed by Fexxv and Fexxvi, since the part of the wind 
seen at low inc hnations is the inn er, highly ionised, surface 
of the outflow. ISim et al. I (l2010bh report emission features 
that are moderately strong, ~140 eV in their example, and 
which are significantly broadened by the Doppler motions 
of the material t o a width of_CT ~ 0.3 keV. This is con- 
sistent with what iNg et al.l (|2010t ) report in a systematic 
study of low-inclination (non-dipping) NS LMXBs with 
XMM-Newton. First, they report an emission line for 80% 
of the sources, i.e. the line is almost ubiquitous. Second, 
they report a weighted average energy for the sample of 
6.67 ±0.02 keV, consistent with Fexxv, and a width of 
(7 = 0.33 ±0.02 keV. Finally they find EWs between 17 
and 189 eV. Inde ed, a qualitati ve comparison o f Fig. 4 of 
Sim et all ()2010al ) with Fig. 3 of iNg et all (|2010l ) indicates 
that even the shape of the lines looks similar. In partic- 
ular, the line of 4U 1705-44 show n in Fig. 3 (third row, 
third column) of Ng et al. ( 2010l ) is strikingly similar to 



the one shown in Fig. 4 (upper-left panel) of ISim et al 
()2010ah . 



As the inclination increases to ^60°, the hottest and 
most ionised layers of the outflow are first intercepted. At 
this inclination the material is almost fully ionised and 
the only important opacity source is Compton scattering. 
Due to the energy dependence of the relativistic Klein- 
Nishina formula for Compton scattering, the opacity will 
be lower and the spectrum will curve upwards at high en- 
ergies. Compared to lower inclinations, the emission fea- 
ture will be slightly broader since a larger colu mn density 



and gra dient in velocity fields will be traversed. ISim et al 



(l2010b h report widths of around 0.35 keV and appear- 
ance of blueshifted Fexxvi narrow absorption lines. At 
even higher inclinations, the part of the wind seen by the 
observer is less ionised and has a larger column density, 
and consequently the absorption features of less ionised 
ions, such as Fexxv, become deeper. The direct compo- 
nent of the spectrum becomes increasingly dominated by 
the scattered/reprocessed component resulting i n a promi- 
nent e mission lines with EWs of up to 400 eV ( Sim et ahl 
2010bl ). This is consistent with the results of this work and 
of previous works on classical dippers (i.e. high inclination 
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sources), in which narrow absorption li nes and broa.d iron 



lines are observed s i multaneously fe.g. ISidoli et al.ll2001 



Parrnar et alJ 2002t Boirin et al 



200911 . 



Boirin et all (120051) and 



■ .20051: iDfaz Trigo et al 
Diaz Trigo et all (l2006ll 



found that a highly ionised absorption plasma is ubiqui- 
tous for these sources and inferred that such an equatorial 
plasma must be present in all LMXBs. Since the plasma 
will develop as a wind or not depending on the size of the 
system and on its luminosity (see Sect. 15 .l]) . in cases where 
the thermal wind cannot develop due to the low luminos- 
ity of the central source the scattering will be produced 
in the hot, static, atmosphere that is revealed through 
absorption lines. 

A caveat of the scenario outlined above may be the 



that iFroning et all (|2008f ) attribute 3/4 of the NIR flux 
to the disc), it is difficult to explain a periodicity in X- 
rays without invoking structure in the accretion disc or 
its atmosphere. Further, the detection of drops of inten- 
sity of up to 80% in the light curves of GX 13-1-1 which are 
energy-dependent strongly resembles the so-called dipping 
behaviour of high-incl ination sources (see e.g. Fig. 1 of 



Dfaz Trigo erahf |2006[) . Therefore, interpretations of the 



spectral/timing behaviour of the source in the context 
of Doppler boosting of a jet pointing directly towards us 
(|Schnerr et al.l 120031 ) are very unlikely at the light of the 
2008 XMM-Newton observations. Instead we favour an in- 
cli nation for the system of ~60-80°. 

Corbet et al. I (^2010^ concluded that the X-ray modu- 



et al. 



relevan ce of radiation pressure in the disc wind 

( 2010allbh considered parameters typical of AGN for the 
disc wind simulations, for which radiation pressure and 
line-driving is co nsidered to be more imp ortant than for 
LMXBs (see e.g. Proga fc KallmanI 2002). For example, 
the inner launching radius is expected to be smaller for 
AGN and the blueshifts of the absorption lines larger. 
However, the large column density of the absorber in 
obs 4 indicates that electron scattering and therefore ra- 
diation pressure is non-negligible and that shielding could 
be even relevant for some parts of the wind. An extrap- 
olation of the reg ion of emission line formation found by 



lation could be caused by dipping behaviour and that the 
structure causing the modulation should not be phase- 
locked. In classical dippers, the dips recur at the orbital 
period of the system and are believed to be caused by pe- 
riodic obscuration of a central X-ray source by structure 
located in the outer regions of a disc resulting from the 
impact o f the accretion flow fro m the companion star into 
the disc (IWhite fc Swanklll982l) . However, the structure is 
locked at phase ^ 0.8 (considering phase as the phase at 
which the disc is eclipsed by the companion). 

We therefore propose a warped, tilted and precess- 
ing disc as the origin for th e varia bility of the X-ray 



Sim et all (l2010b[) of 4.7 x 10^ r^ to GX 13-1-1 yields a ra- modulation. |Dfaz Trigo et al.| ^00^ explained the ap- 
pearance and disappearance of dips in XB 1254—690 as 
well as the changes in the optical emission by invoking 
a precessing warped disc, s imilar to the one in Her X- 
1 ( Gerend fc BovntonI 1976 ). The existence of a warped 



dius of ~10^cm, very close to the estimated launching 
radius of the wind of ~10^° cm. Therefore, dedicated sim- 
ulations for disc winds with parameters appropriate to 
LMXBs are now of utmost importance to allow spectral 
modeling and confirm this scenario. 



5.4. Orbital period dependence of the wind 

We found significant variability both within each observa- 
tion and among all the observations analysed in this work. 
If the disc wind scenario outlined in the previous section is 
correct, we can naturally explain the spectral differences 
as different parts of the wind being probed at different 
phases, i.e. changes in the column density and ionisation 
of the warm absorber/ wind in the line of sight and of 
the reprocessed component drive the spectral variability 
among observations. 

Corbet et al ] (l2010h reported a modulation of ^ 24 d 
based on RXTE/ASM data spanning 14.2 years. They 
also found that the modulation is not constant. For ex- 
ample, their Fig. 6 indicates that for some time intervals 
the dominant peak is 24 d but for others the peak lies 
below (third panel from the bottom) or above (second 
panel from the b ottom) . The fact tha t a modulation is ob- 
served in X-rays (jCorbet et al.ll2010l) indicates most likely 
a high inclination for the source. Although regular dip- 
ping has not been reported for GX13+1, this could be 
due to the long orbital period of ~24 d and the lack of 
a precise ephemeris. While modulation in NIR can be ex- 
plained e.g. by a modulation of the companion (but note 



disc in GX 13-1-1 together with a high inclination for the 
system would explain both the detection of a modulation 
in X-rays and the fact that the structure causing it is not 
phase-locked. Further, the warp would imply that the line 
of sight is slightly different at different phases, explaining 
the fact that we see different parts of the wind along the 
orbital period. The warp could have b een caused by the 
disc wind, as in the case of Her X-1 (jSchan dl' 'l996 l) or 
simp ly by a critical illumination of the disc (jFoulkes et al.l 
20101) . 



6. Conclusions 

We find that: 

- the disc wind in GX 13-1-1 originates at large radii, 
^lO^'^ cm, and is consis tent with be i ng th ermally driven, 
consistent with work by lUeda et al.l (|2004f ) 

- the broad iron line can be explained by reprocessed 
emission in a hot atmosphere, and Compton scattering is 
the main source of its breadth 

- contemporaneous changes of the broad iron line and 
the warm absorber point to an origin in the same disc 
wind. In particular, the broad iron line is too strong to 
be reproduced by refiection at the inner disc if we assume 
illumination from the boundary layer and reprocessing in 
a hot atmosphere, corona or wind illuminated by a harder, 
power-law-like, continuum is needed 
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- the presence of strong energy-dependent obscuration, 
found for first time in GX 13+1, and the absence of eclipses 
combined with the spectral type of the companion, indi- 
cate an inchnation of 60-80°. The emission and absorption 
features expected in the disc wind scenario for such an in- 
chnation are consistent with the observed ones. This gives 
additional support to such scenario. 
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